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Q Abstract. We provide a comprehensive analytic study (rigorous and qualita- 

tive) of the conditions for the existence of a a positive kinematic back reaction 
, term Q > 0, in the context of Buchert's scalar averaging formalism applied to 

^ spherically symmetric Lemaitre-Tolman— Bondi (LTB) dust solutions in which 

^ ^ averaging domains are given as spherical comoving regions containing a symme- 

^_^ try center. We introduce proper volume and quasi-local average functionals and 

f\ functions in order to examine the conditions for Q > 0, and in the process we also 

^yi explore the relation between back reaction, spatial curvature and binding energy 

I for a wide variety of LTB configurations. The back reaction term is positive for 

?-H all "hyperbolic" regular domains with negative spatial curvature, either in the 

OJJ full radial range or in the radial asymptotic range. This result is also valid if 

these domains contain an inner "elliptic" region with positive curvature undergo- 
ing local collapse. For some cases in which positive spatial curvature decreases 
asymptotically, the conditions for a positive back reaction can still be met but 
seem to be more restrictive. Since Q > is a necessary condition for a positive 
^ ,^ "effective" acceleration that would mimic the effect of dark energy (in the context 

\I of Buchert's formalism), we examine this issue in LTB models in a follow up paper 

^~~l (part II). 
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1. Introduction. 

The possibility that observations could be influenced by diferent ways of coarse 
graining and averaging of dust inhomogeneities provides a number of popular 
alternative explanations for cosmic acceleration found in the literature [TJ |51 [31 HJ [51 [B] . 
Among these theoretical proposals, the spatial averaging formalism developed by 
Buchert [Z| and co-workers considers the "effective" acceleration that mimics a 
cosmological constant, which might arise from the so-called "back-reaction" terms 
that emerge by rewriting scalar evolution equations in terms of spatial averages of 
matter-energy density and the expansion scalar. See [Hj for a comprehensive review of 
this formalism and [H] for further discussion. For alternative proposals on averaging 
inhomogeneities, see PUIITT] . 

Unfortunately, whether based on Buchert 's formalism or not, it is technically 
challenging to calculate in practice coarse grained observational parameters or back- 
reaction terms for general non-linear inhomogeneus or "realistic" conditions. Hence, 
besides perturbative approach [T^] and idealized spacetimes [13], proposals that 
examine cosmic acceleration without dark energy have been discussed or tested for the 
spherically symmetric Lemaitre-Tolman-Bondi (LTB) class of dust solutions [U [21 [3 
[m [151 [m [17] • These models [HI [19] are simple, but general enough, inhomogeneous 
spacetimes and so are ideal to test the effects of inhomogeneity. See [501 1211 [HI [2S1 121] 
for further discussion on regularity of LTB models and [19j for a comprehensive review. 

Several articles among the references in [T] [21 [5] [TH [15] have already considered 
the application of Buchert's formalism to LTB models. More recently, Paranjape and 
Singh J^ utilized asymptotic or late time approximations, while Chuang, Gu and 
Hwang [T7] relied on particular cases of exact LTB solutions. These articles apparently 
signal that a positive effective acceleration can occur under certain conditions (negative 
curvature and low density). In the present article and its continuation (part II) we 
enhance and complement this existing literature by looking at this theoretical issue 
in in full general analytic form, without specializing for particular LTB models, and 
without perturbations or approximations. Specifically: we examine in detail sufficient 
conditions for a non-negative back-reaction term, Q > 0, which become necessary 
conditions for the existence of a positive effective acceleration Acs > (applied in 
part II to those cases previously found in this paper to comply with Q > 0). We 
have summarized these results in Table 1 , for the benefit of those readers who wish to 
see them before going into the technical detail. The relation between back-reaction, 
spatial curvature and a binding energy functional, which can be defined in a covariant 
manner [2S], appears in the process of discussing the conditions for Q > 0. This 
relation, summarized in section 14, provides a solid, though model-dependent, support 
for the arguments suggested by Wiltshire [26, on back-reaction, quasi-local energy and 
binding energy. 

We summarize now the section by section contents of the article. We present 
in section 2 the basic parameters and properties of LTB models. In section 3 we 
examine the natural and covariant time slicing afforded by the normal geodesic 4- 
velocity, discussing geometric properties of the hypersurfaces ^T(t) orthogonal to u°, 
which provide the spatial slices in which the proper volume averages associated with 
Buchert's formalism will be defined and computed. Section 4 provides a rigorous 
definition of the proper volume spatial averaging functional, as well as its associated 
local valued function, while in section 5 we introduce the quasi-local average functional 
and function, which will be needed in subsequent sections. The evolution equations 
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for averaged scalars in Buchert's formalism, including the back-reaction term Q and 
its associated effective acceleration A^^g, are given in section 6, while in section 7 
we present the generic sufficient conditions for Q > 0. We discuss the conditions 
for the existence of zeros of radial gradients of scalars ("turning values" or TV's) 
in section 8. The sufficient conditions for Q > are then applied to domains 
in various LTB configurations: "parabolic" (section 9), "hyperbolic" (section 10), 
"elliptic" with open and closed topology (sections 11 and 12) and mixed configurations 
in which an inner elliptic region (enclosing a center) is surrounded by a hyperbolic 
or parabolic exterior (section 13). In section 14 we provide an overview and final 
discussion on the interrelation between back reaction, spatial curvature and binding 
energy. The discussion on properties of LTB models given in sections 2 and 3 is 
complemented with a summary of regularity conditions (Appendix A), the relation 
between radial coordinate and proper radial length (Appendix B) and the fluid ftow 
evolution equations ^7\ UHl UHl EH] (Appendix C) . Appendices D and E provide the 
detailed proof of Propositions 7 and 8. 

The most relevant result in sections 9-13 is the fact that averaging domains 
always exist in which Q > holds for all regular LTB models whose radial asymptotic 
range has, either negative spatial curvature {^TZ < hyperbolic), or radially decaying 
positive curvature (^TZ > elliptic models with open K'^ topology) with sufficiently 
large gradients '^TZ' < 0. In a follow up paper (part II) we explore conditions for which 
Acs > can hold for those configurations in which Q > holds, leading to a robust 
analytic proof on the compatibility between specific LTB models and the existence of 
a positive effective acceleration in the context of Buchert's formalism. 

2. LTB dust spacetimes. 

Lemaitre-Tolman-Bondi (LTB) dust models [HI [19] are the spherically symmetric 
solutions of Einstein's equations characterized by the LTB line element and the energy- 
momentum tensor 

ds'^ = -c^dt^ + ^ dr^ + R^ide^ + sin^ edc/)^), (1) 

^°'' = pc2u°u^ u^^S^, (2) 

where R = R{t,r), R = dR/dr, T = T(r) > and p = p{t,r) is the rest -mass 
energy-density. The field equations reduce to 

R' ^^-f + {T^-l). (3) 

2M' = KpR^R', (4) 

where k = SttG/c^, M = M{r) and R = u'^WaR = dR/d{ct). 

The basic kinematic parameters in LTB models are the expansion scalar, Q = 
Vait", and the spatial trace-less shear tensor aab = Ki^t^ [cUd) — {Q/3)hab 

„ 2R R , , 

^=1^ + 1^' (^) 
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where hat — UaUt, + gab, S"*" = h^'' — Sx^x'' and x° — \/hP' S"^ is the unit vector 
orthogonal to u°- and to the 2-spheres orbits of SO (3). Another important covariant 
object is the electic Weyl tensor E°-^ — UcUdC'^'"^'^ 

E-'^SE^" ^ s^-^£ab^'''^-'^p+^, (7) 

Notice that both aab and Eab can be described by the single scalar functions E, £ in 
a covariant manner. We provide in Appendix C the "fluid flow" evolution equations 
for the covariant scalars p, Q, T, and £. 

The term ^^ — 1 in dsl) is often interpreted as a "binding energy" for comoving 
dust layers [151 HH [53], as ^ is analogous to a Newtonian energy equation. It is 
important to remark that all quantities in pl-Q are are invariant scalars for LTB 
models: M is the conserved "quasi-local" mass of Misner and Sharp [311 [32l l33l l34] . 
R is the "area distance" , while J-" can be related in a covariant manner (in spherical 
symmetry) to a binding energy functional (see section 14 and [2S])- 

It is common usage in the literature to classify LTB models according to a 
"kinematic equivalence class", based on the sign of T^ — 1, which determines the 
existence of a zero of i?^ in (p|. If the given sign holds in the whole regularity domain, 
we have 



J^^- 


-1=0, 


or: 


^=1, 


parabolic models. 


(8a) 


T^- 


-1 >0, 


or: 


J^>1, 


hyperbolic models. 


m 


T^- 


-i<o, 


or: 


- 1 <^< 1, 


elliptic models. 


(8c) 



with the equal sign above holding only at a symmetry center. In general, it is possible 
to consider T^ — 1 changing sign in a given radial range, defining LTB models that 
contain hyperbolic or elliptic domains or "regions" (see [T ^[2^E5] ). 

The standard technique to deal with LTB dust models is to solve the Friedman- 
like field equation (l3| for suitably prescribed functions M(r) and Tir) (the latter 
commonly as a function E = T^ — 1) , using then these solutions to find the remaining 
relevant quantities that may be required for a specific problem. This procedure has 
lead to analytic solutions (mostly implicit or parametric) that are well known and 
have been used abundantly in the literature (see [19' for a comprehensive review) . A 
different approach to study the dynamics of these models is the "fluid flow" or covariant 
"1+3" decomposition of Ehlcrs, Ellis, Bruni, Dunsbury and van EUst PTIESII^IHU] 
(see Appendix C). Standard regularity conditions for these models are discussed in 
[201 mi HH 1131 mi and summarized in Appendices A and B. 

3. Covariant time slicing and geometry of space slices. 

The normal comoving 4-velocity in (fTl) provides a covariant time slicing in which the 
space slices are the hypersurfaces ^T(t\ orthogonal to m", and marked by an arbitrary 
but flxed t = tf). Since we will be considering integral functions and functionals along 
the ■^T(io), it is important to provide a basic discussion on the geometric properties 
of these slices (see Appendices A and B for further discussion) . 

The metric, proper volume element and 3-dimensional Ricci scalar associated 
with the ^T(io) are 

hab^ gab + UaUb^ g^jSl^Sl, ij = r,9,(j>, (9) 

dVp= ^/dct{hab) d^x = J'-^ R^R' sin^ 6* dr d6l d0, (10) 
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mR' R'R i?2 ' ^^^> 

where R = R{to,r), R' = R'{to,r). We will assume henceforth the existence of (at 
least) one regular symmetry center (see Appendix A) marked by r = 0, so that 
R{t, 0) = R{t, 0) = for all t and r > (hypersurfaces '^T homeomorphic to S'^ 
have a second symmetry center at r = re). 

Each ^T(io) is a warped product '^T — x(io,^) Xii:§^(^i0), where the warping 
function is R{to,r) > 0, the fibers are concentric 2-spheres S^ with surface 
area ATrR^{tQ,r), while the leaves x(^0:^) £^re "radial rays" or curves of the form 
[ctQ, r, Oq, ^o], with 9q, (J)q constants. Since the rays at each ^T{t^) are isometric to each 
other, every scalar A function in every ■^T(io) is equivalent to a real valued function 
A : x{to,f) ~^ U^ that corresponds to one of the functions in the one-parameter 
family A{cto,r). Evidently, for time dependent scalars we will have a different radial 
dependence at different '^T(to), while scalars like M{r) or ^{r) have identical radial 
dependence in all ^T(to). 

It is important to mention that the radial coordinate in (IT]) and ([9]) has no inherent 
covariant meaning. However, as we show in Appendix B, as long as standard regularity 
conditions hold, the proper length along radial rays (which are spatial geodesies of 
these metrics) is a monotonically increasing function of r, and so the dependence of 
scalars on r at ant '^T{t) is qualitatively analogous to their dependence on the proper 
length. Since it is evident that all time dependent quantities defined on the ^T{t) 
depend on i as a fixed parameter, we will henceforth omit expressing this dependence 
explicitly. Unless it is needed for clarity, we will use the symbol A(r) instead of A{t^ r). 

The hypersurfaces ^T{t) can be classified in terms of the sign of the spatial scalar 



curvature ''T?. in ( 11 ), and the "kinematic class" given by the sign of T — 1 as in (8a I 



(8cl, which identifies the ^T{t) (or regions of them) as slices of parabolic, hyperbolic 
or elliptic models. Given the existence of (at least) one symmetry center, the topology 
(homeomorphic class) of the ^T{t) are 

• '^T{t) homeomorphic to M? . There is only one symmetry center, at r = 0, hence 



we must have R' > Q everywhere. Notice that this also follows from demanding 
absence of shell crossings (see Appendix A). From the proper length definition 



(2.11 and condition (12. 21) in Appendix B, we must have JF > for these cases. 



which means that this topology is compatible with regions or models of all 



kinematic classes (hyperbolic, parabolic or elliptic). As a consequence of (2.36) 



and (2 Ah), if ^ -^ cxi then i? — > oo, though it is possible to have R -^ constant 
if JF ^ asymptotically. Notice that (2.3&I implies that all ■^T(i) have the same 



topology. 
^T{t) hor 
Since R{t, 0) = R{t, re) = for all t, then there must be a "turning value" r = rtv 



'^T{t) homeomorphic to S"^. There are two symmetry centers, at r = and r — r^ 



of R, so that R'{rtv) — where < rtv < r^- Since the regularity conditions (2.2 1 



and (2.3&I require ^(rtv) = 0, all regular models whose ^T(i) have this topology 
must be elliptic (though elliptic models can also have ^T{t) homeomorphic to 
See section 12 and Appendices A, B and E. 



The interrelation between kinematic class (sign of JF — 1 in (8al-(8cl), radial 



profiles of JF, scalar curvature '^TZ in (11) and topology will be discussed for each case 
in sections 8-13. 
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4. Proper volume average functionals and functions. 

Each regular ||] slice '^T of a parabolic, hyperbolic or elliptic model admitting a 
symmetry center at r = is difFeomorphic to the product manifold 

^ o2//i JN f I^^, M^topology, ,_,„, 

V = ^xS\e,^), ^ = I {, I < , < ,J §3 topology (12) 

Hence, each of these regular ^T contains compact concentric spherical comoving 
regions enclosing a center, bounded by 2-spheres of area 47ri?^(io,?'o)j and 
difFeomorphic to the product manifold 

V[ro]^{r\0<r<ra}cni, (13) 

so that every scalar A in I?[ro] is equivalent to a real valued function A : rjlro] — > M, 
and the set X(I?[ro]) of all scalar functions in I?[ro] is equivalent to the set X(7y[ro]) 
of all real functions in ?7[r'o]. The proper volume of any compact region ^'[t'o] reduces 
to the following real valued integral on ?7[ro] 

"l^p(^o) = / dVp = 47r r' T-'^R^R'dx, (14) 



SO that Vp(0) = 0. 

Definition 1 . The proper volume average functional associated with a 
spherical comoving domain I?[ro] in a given '^T(t) is the linear integral 
functional ( )p : X(77[ro]) -^ M, so that for every scalar function A E X{ri[ro]) 
we get the real number 

Comment . Notice that ( )p is a functional, hence it acts in a non-local manner by asso- 
ciating the real number (A)p[ro] to the scalar A for the whole domain I?[ro]. If we con- 
sider the same scalar A but with its domain given by another comoving region, I'[ri] 
with ri 7^ To, then we obtain by means of ( )p the real number (A)p[ri] ^ {A)p[ro\. 



For all domains ( 13 1 the number tq marks the 2-sphere (fiber) that is the boundary 
of a region I'[ro]. Since rg is a fixed but arbitrary parameter, if it varies then we can 
always construct from the average functional a local valued function of Tq. 

Definition 2 . For every scalar A e X{ri[r]), the proper average function ("p- 
function") is the real valued function Ap : X{ri[r]) -^ M., so that for every 
domain 7/[r] 

Ap{r)^{A)[r], (16) 

where r > encompasses all the domain of regularity of the radial coordinate 



given by ( 12 I 



J We are assuming that the ^T{t) are fully regular, which is not true when a curvature singularity 
arises at a given t. We discuss this issue in Appendix A. 
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Comment . The functional (^)[ro] and the function Ap{r) are closely related, though 
there are subtle but important differences between them. Understanding their simi- 
larities and differences is crucial for the proper understanding of the present article, 
hence we illustrate graphically both objects in figure 1. Notice that for every domain 
r][ro], we have Ap{ro) = {A)[ro], but Ap{r) ^ {A)[ro] for all r ^ vq. 



Notation. In order to avoid confusion and to simplify notation, we will adopt the 
following conventions: we will use the symbol (A) without explicit mention of "[rp]" 
(unless it is needed for clarity) , as it is evident that (A) is the real number associated 
by (151 with a domain I?[ro] defined by (13 1. The functions Ap will be denoted "p- 
functions", reserving the term "average" only for the functional (A). We will use r as 
independent variable of the functions Ap, reserving x or for the integration dummy 
variable. 



Only the functionals (A) can be considered averages of real valued functions A 
taken as a continuous random variables. This follows from the fact that (^4) and Ap 
behave differently under integration along any ri[r]: 



/•r /•?' /-r /-r 

/ {A)Bdx^{A) Bdx but / ApBdx^Ap Bdx. 
Jo Jo Jo Jo 



(17) 



Hence, (A) defined as functionals by (15 1 allow for the construction of momenta such 
as the variance and covariance, as they comply with 



{iA-{A)f) = {A'}-{Af, 
{{A-{A)){B-{B))) = {AB) 



{A){B), 



(18a) 
(186) 
(18c) 



On the other hand, the functions Ap{r) do not satisfy (18a|-(18c), and cannot be 
"averages" of a continuous random variable. The proper volume (14 1 satisfies the 
following properties: 



Vp 
Vp 

v 

V 



(e), 

3i?' {T) 



(19a) 



(196) 



where (19al follows from = [ln(jF ^B?R')\ and the commutation of d/dt with 
the integrals in (15l, while (196l is readily obtained from {F) = 4tt R^ / {3Vp) . The 



commutation rule between the time derivative and the average is 

{Ay-(A) = {eA)-{e){A), 



(20) 



which follows directly by applying d/dt into ( 15 1. The functional (A) and the functions 
Ap comply with the following properties 



{AY 



V 
IP 



[A-{A)\ 



A{r)-{A)[r] = -^ TA'Vpdx, 
'^pV ) Jo 



(21a) 
(216) 



which follow directly by applying d/dr and integrating (15 1 by parts 
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(b) 




x=0 



x=r 



Figure 1. The difference between Ap and (A). The figure displays the 
radial profile of a scalar function A{x) (solid curve) along a regular hypersurface 
•^TYt), together with its dual auxiliary function Ap{x) (dotted curve) defined by 
\1&\ . Panels (a) and (b) respectively display the cases when A' < ("clump") 
and A' > ( "void" ) . The average functional \l^\ assigns the real number (A) [r] 
to the full domain (shaded area) marked by ri\r\ = {x\0 < x < r}, whereas the 
function Ap varies along this domain. Hence, Ap and (A) are only equal at x = r, 
and so they satisfy the same differentiation rules locally, i.e. Ap{r) = {A)'[r] and 
Ap{r) = {A)'[r], but behave differently when integrated along the domain. Notice 
that, as stated in Lemmas 2b and 2c in section 8, and from l |21a| and ( |2f b[ l, if 
A' < in all r)[r] then A — (A) < and the opposite situation occurs if A' > 0. 
This is also true for the quasi-local functions Aq. 
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5. Quasi local average functionals and functions. 

The function M that appears in dsl-Q is for LTB models the Misner-Sharp quasi- 
local mass-energy function, which is an important invariant in spherically symmetric 
spacetimes. It is basically the volume integral of the field equation Q, which will be 
well defined if we assume the existence of a symmetry center (at r — 0) and can be 
given as a proper volume integral "weighed" by T: 

2Af = ^- / pTdVp^K pR^R'dx. (22) 

c Jv Jo 

This integral definition motivates the introduction of a "weighed" average functional 
and its associated local function. These functions will be very useful in the following 
sections. 



Definition 3 . The quasi-local proper volume average functional associated 
with a spherical comoving domain 'D[ro\ is the linear integral functional 
( )q : X{r][ro]) — > M, so that for every scalar A G X{ri[rQ]) we get the real 
number 

CATdVp _ J^^AR^R'dr 
^^^'™ ~ CTdVp ~ J^°R^R'dr ' ^^^^ 

Definition 4 . For every scalar A S X{'D[r]), the quasi-local average function 
( "q-function" ) is the real valued function Aq : X{ri[ro]) -^ M, so that for 
every r € i][r] 

Aq{r) = {A)q[r] (24) 

As with the p-functions, the above definition holds for all 77 [r] and their do- 



main of definition is the full the domain of regularity of r in ( 12 1 



The quasi-local volume averages and q-functions satisfy analogous properties to 



{A) and Ap-. only the functionals {A)q are averages which comply with (18a)-(18cl 
with ( )g replacing ( ). The comparison between (A) and Ap in figure 1 applies also 
to {A)q and Aq. The quasi-local volume is given by 



Vq{r) = / TdVp -- 


= 47r 


Jv[r] 




and its derivatives are 




Vq_m 

Vq- R- ^"^ 




v; 3i?' 




Vq~ R' 





/ R^R!dx^ —R^{r), (25) 

Jo 3 



(26a) 



(266) 



where (26a I follows from O = [ln(_R i?')]' and the commutation of d/dt with the 



integrals in (15) (just as with (|19a[)). The commutation between d/dt and the quasi- 



local average (and q-functions) is exactly the same as (20 1. The following relations 
hold 

A'q^{Aq)'^^[A-Aq], (27a) 

A{r) ~ Aq{r) ^ -^ f A' Vq dx, (27&) 
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which are analogous to (21al~( |21&| . 

Applying the definition ( 23 1 to dsl , (fill and ( 22 1 we obtain the following 
important relations 

K 2M 



Pq 



"R„ = - 



6 [^2 _ 1] 



9^ 



^2 



1 



i?2 



(28 a) 
(286) 
(28 c) 



i?2 

The scalars S and £ in (|6]) and (I?]), associated with the shear and electric Weyl 
tensors, become expressible as deviations of the local scalars 0, p with respect to their 
quasi-local duals: 

S--^(e-e,), (29a) 



^ = -g (p-p?) 



(29&) 



These relations among quasi-local q-functions will be very useful in the forthcoming 
sections. 



6. Buchert's averaging formalism. 

The idea behind Buchert's spatial averaging formalism is to apply the averager proper 



volume functional ( 15 1 to scalar evolution equations for covariant scalars in spacetimes 
under a suitable time splitting (see [3 H]). For the case of LTB models under 
consideration, and given the time slicing furnished by u'^ (see section 3), the functional 



(|15|) must be applied to both sides of the scalar evolution equations: the Raychaudhuri 



and energy balance equations, and to the Hamiltonian constraint (see equations (3.2a 



(3.261 and (3.4| in Appendix C). Using the time derivation rule (20) and (18al-(18c 



we obtain the evolution laws for (0) and (p) and the averaged Hamiltonian constraint 

(6)2 



(6)2 



3 

Poff, 



[Poff + 3 Poff ] : 



9 3 

where the "effective" density and pressure are 

KPeff =k(p) ^ , 

_(37^) 



nPntf = 



Q 

6 2 ' 

and the kinematic "back-reaction" term, Q, is given by 

Q^^((e-(e))^)-6(s2), 



(30) 
(31) 
(32) 

(33) 
(34) 

(35) 



§ These equations follow from the fulfiUnient of the average properties ( |l8a| l-( |l8e[ |, therefore they 
do not involve the p— functions 0p, pp, ^7?.p. 
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where 6S^ = o' nhu "''^ follows from (|6| . The integrability condition between ( 30 1 , ( 32 ) , 
(33 1, (34 1 and (35 1 is the following relation between Q and l^TZ)' 

Q + 2(e)Q + ^(e>f7^) + (='7^>•-o, (36) 

which can be proven to be compatible with the fluid flow evolution equations given in 
Appendix C. 



Equations (30 1 clearly convey the motivation of Buchert's approach: by averaging 



the scalar evolution equations (see Appendix C) we get for a simple dust source 
averaged quantities, (0) and (p), whose evolution mimics that of a source in which 
there is "effective" density and pressure, (33 1 and (34), constructed from "back- 



reaction" terms that arise in the same averaging process. These terms could yield, in 



principle, a form of positive acceleration in (30 1 that arises simply from the averaging 



of inhomogeneities and not from an elusive source like dark energy. 



From ( 33 1 and ( 34 ) , the condition for an "effective" cosmic acceleration that could 



- [pcff + 3Poff] >0, 



mimic dark energy is 

SO that , given a domain of the type ( 
(but not sufficient) condition for (37 

Q>0. 



l{p) >o, 



(37) 



13 1 and p > holding everywhere, the necessary 



IS 



(38) 
Therefore, once we find conditions (necessary or sufficient or both) for the fulfillment 



of ( 38 1 , these will be necessary for (|37[) . We examine these conditions for the remaining 



of the article. 



It is important to notice that both Q and A^jf in ( 37 1 and ( 38 1 can be expressed 
as sign conditions on averages of a single scalar defined in a domain 77 [r] of the form 



Q =3(c*)>o, 

^off = {A.,) > 0, 
with the scalars C* and A given by 

a = C, {x, r) = {Q{x) - (9) [r]f ~ {Q{x) - Q,{x)f . 

2 K 

A = A (a;, r) = -C* (x, r) - -p{x) . 



(39) 
(40) 

(41) 
(42) 



where we used (29a I to express the shear scalar (Tabcr — 6X1 in terms of Q and Q, 



7. Sufficient conditions for a non— negative back reaction. 



Conditions ( 39 ) and ( 40 1 have been examined on LTB models by means of 
approximations J16j or particular exact solutions [15' [T7]. For a more general 
theoretical framework, it is practically evident that finding the general (necessary and 
sufficient) conditions for the fulfillment of these conditions cannot be done without 
resorting to numerical methods, as it requires evaluating average integrals (151 for 



fully general metric functions, like R and R' , which are known (at best) in implicit or 
parametric form from solving (pi. 



However, if we are interested in finding sufficient conditions for (37l-(40l, it is 



(fortunately) not necessary to compute integrals like (15). If what is needed is simply 
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to find out the sign of an averaged quantity (A) in a given domain 77 [r], we can obtain 
sufficient information on this sign simply by looking at the behavior of A point by 
point in the domain. Concretely, we will use the following property of integrable 



functions A defined in domains like (13 1: 



A{x)>0 V xeT][r] ^ {A)[r]>0, (43) 

though, it is important to mention that the converse is not true: we can have (A) [r] > 
even if A < holds in subsets of r][r]. We will use the property (43 1 to examine 
sufficient conditions for ( 39 1 , and then look at their implications for ( 40 1 . 



Sufficient (but not necessary) conditions for ( 39 1 and ( 40 1 are given by 
C*(a;,r)>0, 
A*{x,r) > 0, 



(44) 
(45) 



holding for all x g r][r]. Notice that if we prove (44 1 and/or (45 1, then Q > and/or 
^eff > follow automatically from ( 43 1 as a corollary. 



Testing conditions (44|-(45) can be very difficult because C^,{x, r) is really a fam- 



ily of functions of x for arbitrary fixed r. Fortunately, these conditions can be greatly 
simplified by means of the following: 



Lemma 1 : (W) = in every domain ri[r] for W — W{x, r) given by 

W{x, r) = [Q{x) - (6) [r]f - [Oix) - Opix)]' , (46) 

Proof . Expanding ( [46| and applying ( [T5| we obtain with the help of (18&I 

(W)[r] = -(e)[r]2 + ^ /' [266^ - 6^] V'^dx. (47) 

•^pV) Jo 

Inserting 6 = V'p/Vp and 6p — Vp/Vp in the integrand above, and bearing 
in mind that (0) and Op coincide at the domain boundary x — r, leads to 
the desired result: 



(>V)[r] = -(e)[r] 



Vp(r) Jo 



dx 



-(e)[rp + e^(r) = o. 



(48) 



An analogous result follows for the quasi-local average acting on a scalar like 



W with ( )q and Qg instead of 



and Qr. 



Corolary: for any domain ri[r] we have Q = {C^)[r] = {C)[r], with C = C{x) 
given by 

c = [6 - epf - [6 - e,]' , (49) 

The proof follows directly from Lemma 1, as {{Q — (0)[r])^) ~ ((O — 0p)^). 



By using Lemma 1, sufficient conditions for Q > given by (39) can be rewritten now 



in terms of C, which can be given as a function of r (since x is a dummy variable): 

Cir) - [e - Opf - [6 - &gf 

= [6, - Gp] [6 - e, + e - Gp] > ^ Q> 0. (50) 
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The behavior (sign) of this quantity must be examined for every domain 77 [r]. 
Considering (216) and {27 b\ applied to Q we get 

e{r) - Bpir) = -^ [ e'{x)Vp{x)dx, (51a) 



e(r)-e,(r) 



Vp(r) Jo 
1 



V,ir) 



e'{x)Vq{x)dx, 



(51&) 



and inserting these expressions into (50 1 and rearranging terms, we can express this 
condition as 



C{r) = $(r)*(r) > 

$(r) = f e'{x)Lp{x,r) 
Jo 



Q>0, 



dx. 



*(r) = 
with if and ^ given by 

ip{x,r) = 

■ip{x,r) = 



Q'{x) ip{x,r) dx, 



Vp{x) Vq{x) ^ Vp{x) 

Vp(r) V,(r) Vp{r) 

Vp{x) Vq{x) ^ Vp{x) 

Vp(r) V,(r) Vp(r) 



Tp{x) 



J^pjx) 



where Tp is the p-function associated to T , and we have used the relation 



V,(^) 
Vp{x) 



J~pyX)^ 



Vp(r) 



•^p(^) 



(52) 
(53) 

(54) 

(55) 
(56) 

(57) 



which follows directly from (15 1 and (23). Since we need to compare the values oi Tp 
in interior points x with those in the boundary of 77 [r] (for arbitrary 77 [r]), the following 
expressions will be very handy: 






y,T-Tp 



Va 



T 



3R J^ J^r, 



R 



T{x) ~ Tp{x) 



T'{x)Vp{x)dx, 



(58a) 
(58&) 



'^p{^) Jo 

We remark that these properties are valid for every ^T(i) and for all x G r/lr]. Notice 
that, even if J-" = 0, the function J^p involve Vp and Vq and so is time dependent. 
However, the radial profiles of J-p in all ^T{t) are qualitative analogous (see Appendix 
B), while the profile of JFp is analogous to that of JF (see Lemmas 2a-2c in the following 
section) . 

The fulfillment of ( [44| ) is now equivalent to that of ( [52| ), and it clearly depends 
on the signs of ip and "0 (besides the sign of O') at all points in any domain along 
arbitrary ^T(i). Since, by their definition, Vp(0) = V,(0) = and T{0) = 1, ([56]) and 
( 55 1 imply 

^-(0,0 = 0, 7A(r,r) = 2, 



d_ R^{x)R'{x) Tp{r) + Tix) 

gj{x,r)- ^^^^ Vp{r)Tp{r) ' 

V5(0,r) = <^{r,r) = 0, 
d i?2(x) R'{x) Tpir) - T{x) 

dx'^^'"'''' ~ :F{x) Vp{r)Tp{r) ' 



(59a) 
(596) 

(60a) 
(606) 
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which indicates that as long as J- and i?' arc non-negative, the signs of ijj and dip /dx 
are non-negative, and so ^ basically depends on the sign of Q' . On the other hand, 
the signs of Lp and dip/dx are not determined, hence the sign of <1> requires more 
examination as it depends on both: the sign of & and the ratio !Fp{x) / !Fp{r) (which 
will depend as well in some cases on the sign of T') . In elliptic configurations where 
the ^T(i) have spherical topology, !F and R' can become negative, and so the sign of 
^ is also indetermined (see section 12). Since the sign of C depends on the sign of the 
product <i>^, we will need to obtain the conditions for both terms having the same 
sign. 

8. Radial profiles and turning values ("TV's"). 



The fulfillment of conditions (44) and (45 1 (or ([52|) is strongly dependent on the radial 



and Vq. We will 



profiles of scalars, such as i?, T, Tp, 6, 6p, Qq and the volumes Vp 
need to probe these profiles along domains rj[r] in the space slices ^T of various LTB 
configurations. In particular, it is important to examine the cases in which radial 
gradients can vanish in a domain 77 [r]. 



8.1. Turning values (TV's). 

Definition 5 . We will denote by "Turning Value of a scalar A" (TV of A) a 
value X = rtv G ri[r] such that A'(rtv) = under regular conditions (which 
excludes shell crossing singularities and surface shells, see Appendix A). 

It is important to remark that the TV's for each of i?, T and Q occur under different 
conditions. In general, a given rj[r] could exhibit either (or all of) these TV's in 
different coordinate values of r. Profiles and TVs of r-dependent functions like J-{r) 
and M{r) can be directly associated with initial conditions at a fiducial ^Ti = '^T{ti) 
(see Appendix B), hence they are marked by the same comoving radial coordinate in 
all ^T{t). TV's of R only occur in elliptic models in which the '^T(i) have topology 
§^, and so will also be marked by the same value of r. In this case, as shown by 
the regularity condition (2.2 1, the TV must coincide with a zero of !F (see section 



12 and Appendix B). On the other hand, the TVs of other scalars (O, p, TZ) are 
not (necessarily) related to initial conditions, hence their coordinate location will (in 
general) change from one ^Tit) to the other. The exception to this rule occurs when 
there is a TV of R, as in this case the TV is common to all other scalars save T (see 
section 12 and Appendices A, B and E and for further discussion). 



8.2. TV's of T and spatial curvature. 

The quantity 1 — T^ and its gradients are directly related to the spatial curvature 
^TZ in (111 and ca n be associated with initial conditions. It is useful to examine its 
relation to ^TZq in (28&), its quasi-local function obtained from (23 1. From (111, (23 1 



and (28&) we obtain 



T' 



\"RqR- 



T' 



RR! 



"R 



_ 3 

'7^„ 



^TZR^R'dx, 



RR' 



'Tla 



2R'/R 



(61a) 



(61&) 
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where we used (|27a|) and (|27&| to eliminate '^TZ' in terms of R' 



^n and '^Uq. Equa- 



tion ( |61&P will be used when discussing elliptic models or domains in which TV's of T 
occur (though notice that T and T' must also comply with the stringent constraints 



given by (1.3 1 and (1.4a|-(1.4c)) 



8.3. TV's of e. 



The condition for such a TV (with R' > 0) come directly from the constraint (3.3 1 



(Appendix B), which with the help of (29a I can be rewriten as 

'3R' 



e'^Qo 



R 



.^(0)(l + j(e)) + (^5(e)y 






eye, 

3R'/R ' 



(62 a) 



{62b) 



Since 0„ is related through (28c) to M, T and R, it can be expressed also in terms of 



(second order) radial gradients of these functions. A similar condition can be obtained 
for the TV of p through (p^l) and (p9a|-(p96l). Examples of LTB models with TV's 



of p and 8 are found in [35, ^ 



S.J^. TV's ofAp and Ag. 

We examine in the following Lemmas how the radial profiles and TVs of scalars A^ 
and Aq relate to the radial profiles and TV's of A: 

Lemma 2a . If there is no TV of A in ri[r] then there are no TV's of Ap nor 

Aq. 

Lemma 2b . If A' > for all x e ri[r], then A{r) > Ap{r) and A{r) > Aq{r). 
Lemma 2c . If ^' < for all x e r][r], then A{r) < Ap{r) and A{r) < Aq{r). 



Proof. These results follow directly from (21 6 1 and (276 1. The converse state- 



ments are not true (for 77 [r]). Notice that these results are valid for integration 
domains like ri[r] in (13 1, which contain a symmetry center. Figure IT] illus- 
trates these Lemmas. 



Lemma 3 . If there is a TV of A at x = rtv and i?' > in 77 [r], then for 
sufficiently large r there will be a TV of Ap a,t x = ri > rtv and a TV of Aq 
a,t X = r2 > rtv, with A{ri) = Ap{ri) and A{r2) = Aq{r2). 

Proof. Let A' pass from positive to negative at a; = rtv • As x reaches rtv the 



integral in (21&I is still positive and so A(rtv) > ^p(rtv), but for rtv < a; < r 



the integrand becomes negative, and so the contributions to the integral are 

increasingly negative. Since V{x)/V{r) is increasing, if r is sufficiently large, 

then a value x = ri > rt 

vanishes (thus A{ri) 

for X > ri 

positive. The proof is identical for Aq, but using (|276| instead of (|21&l 



is necessarily reached so that the integral in (21&I 

Ap{ri)). From ( |21a| , we have A^(ri) = and A'^ < 

An analogous situation occurs when A' passes from negative to 

If 




i?' = in ri[r] (there is a TV of R), then we can have A' = and A' 
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with Ap ^ A and Aq ^ A. The result of this Lemma is illustrated in figures 
HI It] and [Hj where we compare the profiles of T and !Fp . 



9. Parabolic domains and models. 



From (8a I: J-" = 1, J-"' = for all 7/[r]. Hence, from (111 we have TZ = at all a; G ri[r] 
and so these domains are spatially fiat. The topological class is necessarily M^. 

Proposition 1: C = holds in all regular parabolic domains 77 [r]. 



Proof. 

If JF = 1 holds for all x £ ?y[r], equations (58al and (586l imply Tp = 1. Hence, (p{x,r) 

in (55 1 vanishes identically in ri[r], and so <i> and C also vanish. Likewise, ii J^ = 1 for 



X e T][r], then Ap — Aq and so C is identically zero. 



Note. 

It is important to remark that the converse of Proposition 1 is false: C — Q holding 
in a given domain •q[r] does not imply that the domain is parabolic. It is possible to 
obtain a vanishing ip for specific cases with T ^1. This possibility was overlooked in 

[Tsmz]. 



10. Hyperbolic domains and models. 



From (8&I, hyperbolic domains and models are characterized by J^ > 1, with !F{Q) = 1. 
The topological class is necessarily R'^. From (2.2) and also from demanding absence 



of shell crossings [20j |2T1 |22j |23l [24] (see Appendix A) we must necessarily have 

T'>Q VxeryH, (63) 

with T' — Q only at r = (or possibly as r — » 00). Regularity conditions require also 
absence of TV's of R. From (11) and (2.2) we have 

3i 



hyperbolic and regular domain 



^7^<0 Va;e77[r], 



(64) 



but the converse is not true, as '^TZ can be negative in regions of elliptic models (see 
section 13 and figure ]9l. However, from (286) and (61al-(61&l, for the quasi-local 



curvature scalar ^TZ-q we have 



hyperbolic and regular domain <4> iZq < Va; € ?7[r]. (65) 

Regularity conditions [22], ^23^ do not allow for a hyperbolic region (enclosing a center) 
to be surrounded by parabolic or elliptic "exteriors" . Thus, as long as regularity 
conditions hold, the existence of regular hyperbolic domains ry[r] implies a full regular 
hyperbolic model, with full asymptotic radial range, in which '^TZ < 0, .?^ > 1 and 
.F' > hold along all the ^T{t). Since regularity conditions prevent TV's of R and T, 
the only possible TV is that of Q. We examine first the the case without TV's and 
then the case of a TV of O 
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10.1. Hyperbolic domains without TV's. 

Proposition 2: C > holds in all regular hyperbolic domains r][r] without a 

TV of e. 



Proof . 

From ( 63 1 and Lemma 2a (previous section) : 

T'>0 Vx e 7][rl 



(66) 



hence for every r][r] we have Tp{x) < Tp{r), which together with (58al-(58&l and 
(59al-(60&) implies that ip and ip in (551 and (56 1 are both non-negative (see top 
panel of figure [2]). Then, if & > 0, we have $ > and \E' > 0, and their product 
is non-negative. If O' < 0, then $ < and vp < and their product C is non-negative. 



Note . 

Proposition 2 holds for all domains ri[r]. In this case, C > (and consequently) Q > 

hold for all regular hyperbolic models with monotonical profile of Q. 



10.2. Hyperbolic domains with a TV of Q. 



If there is a TV of Q (see (62a|-(62&l), then the integrands of $ and ^ no longer have 



a determined sign. The fulfillment of C > becomes domain dependent. Considering 
(58al-(58&) and (59al-( [606| ), this situation is illustrated by figurel2] whose lower panel 
shows the profiles of Q' {x)ip{x,r) and Q' {x)'4){x,r) when 0' changes sign in a given 
ri[r]. We consider this case in the following 

Proposition 3: In regular hyperbolic models with a TV of Q there always 
exist domains 77 [r], with < r < ri or r > r2 with ri < r2, for which C > 
holds. 

Proof . 

Let X = Ttv mark the TV of 0, so that r/rtv > 1 (all domains with r < rtv have 
no TV's and so are equivalent to those of Proposition 1). We consider the case in 
which 0' > for < x < rtv and O' < for r^^ < x < r (see bottom panel of 
figure [2]). The case with opposite signs is analogous. Notice that each of the functions 
ip{x, r) and ip{x, r) in ([55l)-([56| is really family of functions of x with r appearing as a 
different fixed parameter for each domain 77 [r]. As a consequence, the curves shown in 
the bottom panel of figure [2] (and so the magnitudes of {a, /3, 7, S}) will be different 
for each domain rj[r] in any ^T. However, since these curves must all comply with 



(59al-(606) in all domains, they have the same qualitative form as those shown in 
figure [2|in all T^i?"]. [jj] 

For curves as those in the bottom panel of figure [2] we can write ( 53 1 and ([54| as 



/•rtv rr rrtv rr 

$= / eVda;+ / e'ifdx^ / B'lpdx- / \e'\ipdx, 

Jo Jrtv Jo Jrtv 

/•rtv rr rrtv rr 

*= / e'V'da;+ / e'^|Jdx^ / e'^pdx- / |e'|V'da;. 

Jo Jrtv Jo Jrtv 



(67a) 
(676) 



II This qualitative equivalence of the profiles of ip{x,r) and ip{x,r) in all r][r] (with or without the 
existence of TV's) applies to the proofs of all subsequent Propositions. 
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x = 




x = 



x = rt 



x = r 



Figure 2. Functions ip and i/> for hyperbolic domains. As shown by the 
top panel, both functions ip{x,r) and ip{x,r) are non-negative, hence the signs of 
$ and ^ depend on the sign of ©'. The lower panel displays the functions ©' (/p 
and 0' ip (integrands of <I> and >!') when there is a TV of and 0' passes from 
positive to negative at x = y inside ri[r]. For each profile of the integrands we 
can select four numbers {a, /3, 7, S} to constrain the integrals. As we show in 
Proposition 3, it is always possible to find domains »)[r] for which C > 0. 



For all such integrals, we can always find four positive real numbers {a, (3, 7, 5} such 
that (see bottom panel of figure l2|: 



< / eVda;< artv, < 



0< 



Q'ipdx < 7rtv, 



0< 



|eVdx</?(r-rtv), (68a) 
|e'|V'dx<(5(r-rtv). (686) 



Given (67aH(67&l and (68aH(68&l, it is straightforward to show that C == $* > 



holds for all domains rj[r\ such that 



< min 



1, 



7 
5 



1 



(69) 



07 r a 

hold, while C — <i>^ < holds for intermediary complementary values of r/r^^. The 
inequalities ( 69 ) provide the required values of ri and 7^2 for the domains for any given 



profile of OV and OV for a single TV of 0. 



Notice that the value of rtv follows from the fulfillment of 0' = in (62 a I and is 



fixed at each '^T, irrespective of the domain 77 [r] that we might choose. This means 
that in domains with r > rtv, but r close to rtv, the contribution of the areas cor- 
responding to < a; < rtv will dominate the contribution from rtv < x < r (which 
accounts for inner zone of the fulfillment of C > 0) . Since we can choose a domain 
as extended as we wish [i.e. r can be as large as we want), the curves of Q'lp and 
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Q'lp for each increasing domain will be like those in figure p] but with the ratio rt^/r 
becoming smaller as r grows. Since there is an asymptotic range of r, then for a suf- 
ficiently large r the contribution of the area corresponding to rtv < x < r will easily 
overtake that corresponding to < x < rtv (which proves the fulfillment of C > in 
the outer asymptotic area). Hence, domains always exist so that C > holds in in the 
inner and outer ranges given by the Proposition. 

Note . 

The most important implication of Proposition 3 is the fact that a TV of & merely 
leads to the existence of an "intermediary" zone with negative C, but no matter how 
large can this zone be, we can always have C > in the appropriate ranges. If @ has 
several TV's, then there could be a complicated pattern of intermediate zones, but 
C > will always hold (at least) in the asymptotic range as long as there is a clear 
asymptotic sign of Q'. 

11. Elliptic domains and models with topology M'^. 

The regularity condition (12. 21) implies absence of TV's of R, so that we have a full 



asymptotical range of r. The range of JF in (8cl reduces to < .7-" < 1, with J-{0) = 1 
and JF = only (possibly) as r ^ oo. There are either no TV's, or TV's of O or 
J^ or both (only restricted by regularity conditions). Spatial curvature is positive 
everywhere as long as there are no TV's of JF. We examine first the case without TV's 
and then cases with TV's. 

11.1. Domains without TV's. 

If there are no TV's of !F and 0, then 

T'<0 VxeryH, (70) 



with J-' — only at r = (or possibly as r — > oo). Equations (61al-(61&) imply that 



holds everywhere in ?/[r]. Elliptic domains without TV can arise in the following 
situations: 



• 



Near a symmetry center (see Appendix A). Condition (71 1 will hold for all elliptic 
domains ri[r] with r w 0, since ^U w ^TZ^ w ^7^(0), 6' w and T' < 0. 



Curvature voids. From (61al-(61&) and Lemma 2b, condition (71) holds if we 
have ^TZ' > for all x £r][r]. 

Near homogeneous curvature. If curvature gradients ^TZ' , ^TZ' are not large 



enough, then equations (27aH(27&) imply "^TZ ~ TZ^ and condition (71) holds 



3-0 ~ 31 



However, even in arbitrary elliptic regions or models with turning values of J-", & and 
R, marked by (say) ri, r2, r^, all domains r][r] with r < min(ri, r2, r^) will be domains 
without TV's. 

Proposition 4: C < holds in all regular elliptic domains ri[r] without a TV. 
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From ( 70 1 and Lemma 2a, we have 
< J^„ < 1 



t;<o. 



20 



(72) 



Therefore, for every ri[r] we have J-p{x) > ^p{r), which together with (58a)-(58&l and 
(|59a|-(l60l implies that </? < and V^ > 0. Then, if 6' > 0, we have $ < and * > 0, 



and their product is non-positive. If O' < 0, then $ > and ^ < and their product 
C is non-positive. 

Note. 



Proposition 4 shows that eUiptic domains without TV's and complying with ( 71 1 



have negative back-reaction and so they necessarily exhibit an effective deceleration. 



Provided the regularity constraint (1.3) (Appendix A) is satisfied, this result becomes 
domain independent in a given '^T(t) if there are no TV's for all domains r][r] in it. 
However, it is extremely likely that in any regular elliptic model the ^T(i) will be free 
from TVs only for a reduced range of t. 



x = 





x = 



Figure 3. Functions (p and tp for elliptic domains without TV's. The top 

panel shows ip{x, r) < and <p{x, r) > for this type of domains. The lower panel 
displays the functions ©' if and 0' ip (integrands of <I> and 'I/) when ©' passes from 
positive to negative at x = y inside rj[r]. 



11.2. Elliptic domains with a TVoJQ. 

If there is a TV of Q (but no TV's of T and R), the signs of <& and ^ are no longer 
determined and the fulfillment of C > becomes domain dependent. Considering 



(58al-(58&) and (59a|-(606l, this situation is illustrated by the bottom panel of figure 
|3J displaying the profiles of Q' {x)Lp{x,r) and Q'{x)'4!{x,r) when Q' changes sign in a 
given ?y[r]. We consider this case in the following 
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x = Q 



x = y 



x = r 



Figure 4. Elliptic domains with a TV of J^. The figure displays the function 
y^{x) (thick curve) when it has a TV, together with its dual function Tp(x) (thin 
solid curve) in a domain ri{r\ (shaded region) in a given hypersurface ^T{t). The 
relation between the profiles and TV's of J^{x) and Tp{x) is given by Lemma 3. 
Notice that the relation between Tp{x) and Tp{r) depends on the value x = y, 
as shown in ([75b. The location of y is (in general) different from that of the TV's 
of JP" and J^p and it depends on the choice of r. From l |55[ l, the value x = y 
divides every ri[r] in two zones; one {0 < x < y) with ip{x, r) < and the other 
(y < X < r) with ip(x, r) > 0. 



Proposition 5: If there is a TV of in a regular elliptic region at x = rtv, 
without TV's of J-' or R (open topology), then domains 77 [r] exist for which 
C > holds for ri < r/rtv < ?'2 with ri, r2 > 1. 

Proof . 

We consider the same sign configuration for Q' as in Proposition 3. We use (67a I 



(67&I and (68a)-(68&l to constrain <!> and ^, bearing in mind that now ip < (see 
figure Isl. We readily obtain that C = $^ > holds for all domains r][r] such that 



a 



^,l,^ + l|<-<max|- + l,- + l|, (73) 

Given a profile of (with one TV) in a domain ri[r] and its associated ratios a/P 
and j/S, condition (73) will always hold for elliptic domains in hypersurfaces '^T{t) 
allowing for an asymptotic range (topology R^), where the ratio r/rtv > 1 can be as 
large or small as required. Notice that Proposition 5 is not affected by the fact that 
curves in figure^ are different for each 7y[r], as they are qualitatively analogous and r 
is arbitrary (see comment after Proposition 3). The case with the ^T{t) having closed 
topology is discussed in Proposition 8 and figure [6] 



Note . 

Since elliptic models allow for a collapse of initially expanding dust layers, in general, 
Q for any fixed r will pass from 00 (initial singularity) to — cx) (collapsing singularity), 
with a likely complex radial dependence pattern along the '^T{t) at all times. Thus, 
even if we assume that there are no TV's of JF or R, it is extremely likely that the 
^T(i) will exhibit at least a TV of G. Proposition 5 shows that such a TV introduces 
for some elliptic models with open topology (which satisfy (71l) an "intermediary" 
range in r with positive back-reaction, while back-reaction remains negative in the 
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region around the center and also in the asymptotic range. As we show below, if there 
is a TV of J-, then back-reaction can be positive in the full asymptotic region. 



11.3. Elliptic domains with TV's oj T (but not Q). 



We consider first the case when there is only a TV of JF. As shown by (61al-(61&), 
the necessary condition for a TV of JF is "^TZ'.^, < 0, which (by virtue of (27al-(276)), 
implies a negative gradient ^TZ' . If the gradients "^TZ' and "^7?.'^ are sufficiently large 
and negative (assuming regularity), then ^TZ and '^TZ^, will be sufficiently different from 



'7^(rt,) = 



(74) 



each other (from (27a|-( p7&| )) to allow for a "critical value" 

'7^,(rtv) 
3 ' 

to occur, where rtv marks the TV of T. Hence, we have JF' < (or '^7?. > '^TZ^,/3) for 
< r < rtv and J^' > (or ^7Z < ^TZ^/3) for < r < rtv (see figure |4|. 

While a function !F with a TV can be prescribed as an initial condition, it is 
important to remark that J- and its gradients T' cannot be arbitrary, as they are 



strongly constrained by regularity conditions (see equations (1.3l and (1.4aH(1.4cl in 
Appendix A). 

Proposition 6: If there is a TV of JF in a regular elliptic region, but without 
a TV of 0, then domains 77 [r] exist for which C > holds for r > a for some 
a >0. 



Proof . 

Considering Lemma 3 and (58al-( [58&| , the existence of a TV of J^ implies a profile 
of J^p of the form illustrated by figure [4] Bearing in mind (551, there always exist 
y G ri[r] (in general ^ rtv) in every r][r] and in every ^T{t) such that 



Tp{x) > Tp{r), 
Tp{x) <Tp{r), 






Q<x<y, 
y < X < r. 



(75) 

The TV of T has no effect on the sign of V', so the profiles of Lp and ij) take the forms 
shown in the top panel of figure [5] The sign of C then depends only on the sign of 
$, whose integrand ip changes sign due to the presence of the TV of J- (in agreement 



with (751) 



We use (58a|-(586), (59a)-(60&l and (75) to set up an analogous construction as 
in (67a) and (68al applied to ip and ip displayed in the top panel of figure p] Since 
^ > and there are no TV's of O and i?, we only need to deal with <I>. Considering 
the numbers a, /3 to constrain <i> for the sign change oi ip at x = y (see figure HI), it is 
straightforward to show that C > holds for all domains 77 [r] such that 



r a 

y-^ 



1, 



(76) 



where the ratio a/P depends on the profile of ip(x,r) and on the value y such that 
(p{y,r) — 0. As in Propositions 2 and 5, the profiles of Q'lp and Q'lp (as displayed in 
ffl) will be qualitatively analogous in each domain 77 [r], hence it is always possible to 
satisfy a condition like (76) because there is a full asymptotic radial range and so r/y 



can be as large as needed. 
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11. 4-- Elliptic domains with TV's of T and Q. 

The signs of $ and ^ now also depend on the sign of 0' and we have the following 

Proposition 7: If there are TV's of T and O in a regular elliptic region, 
but without a TV of i?, then domains T][r] exist for which C > holds for 
< ri < r and r > r2 with ri < r2. 

Proof . 

If there is a TV of marked by a; = rtv (besides the TV of T) , the effect is simply to 
make ijj change sign once (at the TV of Q), with an extra sign change in Q'ip. This is 
illustrated by the lower panel of figure [5J The detail of the demonstration is given in 
Appendix D, which provides ri and r2 in terms of y/rtv and r/r^v for any given pro- 
files of 0'(y9 and Q'ip. As in previous Propositions, provided the regularity constraint 
(1.3 1 (Appendix A) is satisfied, the fact that the '^T(t) have open topology allows for 



ratios r/rtv and y/rtv as large or small as needed to fulfill C > (see Appendix D). 

Note. 

The effect of TV's of T and Q is to allow for positive back-reaction in the asymptotic 
range away from the center, in elliptic domains in which '^TZ shows a sufficiently large 
decreasing gradient and when hypersurfaces ^T{t) have topology M'^. We emphasize 
again that !F and !F' must comply with regularity conditions (see Appendix A). 

12. Elliptic models with topology S'^. 

In this case the radial range is restricted by: < x < re, where x — re marks the 
second symmetry center. Since R{t,0) = R(t,rc) — 0, then there must exist a TV 



oi R dX X — rtv. The regularity conditions (2.2 1 and (2.3&I necessarily imply that 
^{rtv) — 0, and so 

R'>0, 0<J^<1, 0<a;<rtv, 

R' <0, -1 < J^ < 0, rtv < a; < ^v, 

(77) 

so that !F' < for all x E ri[r], with T'{0) = ^'{rc) — 0. For elliptic models of this 
type, this TV of R is also a TV of the scalars Q, p, ^TZ, M and all Aq, but not T and 
J^p (the proof is given in full in Appendix E). However, in locations r^ y^ rtv there 
could exist other TV's of !F (as long as T{rtv) — holds) or 6. Depending on the 
combination of these TV's, some domains rj[r] can be analogous to the domains of the 
type examined in Propositions 4, 5, 6 and 7. 

Lemma 4: In any ^T{t) of an elliptic model with topology S"^ in which 
R' {I'tv) — J^i^tv) = 0, there always exist a value x — y > rtv so that every 
domain r][r] with r < y is analogous to a domain of a '^T{t) with topology 
M^' . We have then the following types of domains for 77 [r] with r < y: 

• If y < rtv and there are no TV's of JF or O in 77 [y] , then ri[y] is equivalent 
to an elliptic domain without TV's, as those considered in Proposition 
4. In particular, all domains sufficiently close to the center will be of 
this type. 

• If 2/ > rtv and there are is no TV's of !F, then there is only the common 
TV of R and at a; = rtv In this case, ri[y] is equivalent to an elliptic 
domain with a TV of Q, as those considered in Proposition 5. 
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x = 



x^y 



x = r 



x-rt 



&>0 


0'<O 


3^ ^ 




X 


&^ ^^ i 


"< 




(p<0 ; 


cp>0 



x^y 



X = r 



Figure 5. Profiles of integrands of <J> and 'I' in elliptic domains writh 
TV's of ^ and 0. The top panel displays (j> and i/i defined by l |55| and ( |56[ l 
when there is a TV of J^, while the lower panel shows @'4> and Q'lp when there 
are TV's of JF and (at r'tv) with 0' passing from positive to negative. 
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Figure 6. Elliptic domains with a TV of R. The figure displays the function 
J^{x) (thick curve) together with its dual function J^p{x) (thin solid curve) in a 
domain r][r] (shaded region) with a TV of R marked by a; = rtv (which is then 
also a zero of J^ and a TV of 0). Notice that J'-" < for r > rtv and J^p < 
for r > y. Hypersurfaces ^T{t) have S^ topology with a second symmetry center 
at Tc- Notice that all domains with r < y are analogous to those examined in 
Propositions 4, 5, 6 and 7, whereas in Proposition 8 we considered domains with 
r > y. 
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• If y < Ttv and there is a TV of JF at z 7^ rtv but no TV of 0, then 77 [y] 
is equivalent to an eUiptic domain with a TV of JF, as those considered 
in Proposition 6. 

• If y > rtv and there is a TV of J^ at z G ri[y], then 77 [y] is equivalent 
to an elliptic domain with TV's of J^ (at a; = z) and (at x = rtv), as 
those considered in Proposition 7. 

Proof: Since JF is monotonous, then Lemmas 2a and 2c imply 



•^;<o 



and J-p > J- Vx e 77 [r] 



(78) 

Thus, there will always exist y G 77[r] such that y > rtv and ^p{y) = hold, 
so that Tp{x) > for all x < y. The profiles of JF and !Fp are shown in figure 
|6] These profiles and the fact that x ~ rtv necessarily marks a TV of Q 
imply that the above mentioned combinations of TV's are always possible. 



Note: While Propositions 4-7 apply to all these situations, the restricted 



radial range < a; < y implies that conditions of the form ( 73 H( 76 ) or (4.3 1 



(4.4 1 could fail to hold (not to mention failure to comply with the regularity 



constraint (1.3), see Appendix A). Hence, when these domains are part of 



an elliptic model with closed topology, the wording of Propositions 4-6 must 
be modified by stating that "regular domains complying with C > might 
exist" . 

If we consider domains ri[r] with r > y > rtv, where ^p{y) = 0, we have a situation 
without analogue in elliptic models with open topology because J^ and J-'p are both 
negative (see figure [6]). Assuming that the only TV is the one marked by rtv (that is 
a common TV of 0), we have the following 

Proposition 8: If there is a TV of i? at a; = rtv in an elliptic model, without 
other TV's in r 7^ rtv but with J-p{y) = for y > rtv, then domains ri[r] with 
y < r < re might exist for which C > holds. 

Proof . J-' and J-'p monotonously decrease passing from positive to negative (as 
displayed in figure l6]). Since T < Tp but —1 < Tpir) < 0, Lemma 4 implies that 
there is always a value y > rtv such that J-p(x) < for y < x < r. The functions 



ip{x,r) and tp{x,r) in (55 1 and (56 1 take the forms: 
Vp{x) 



V 



V 



Vp(r) 
Vp(x) 



1- 



J^pjx) 

l-^pWI 



^ 



^ 



Vp{r) 
Vpjx) 
Vp{r) 



1 



\Mr)\. 
\^p{x)\ 
\^pir)\ 



0<x<y, 



y < X < r, 



(79 a) 



(796) 



Since |.Fp(r)| < 1 and |^p(a;)| < |.Fp(r)| in y < a; < r, then (p{x,r) > in all ?7[r]. 
On the other hand, there is also a value z < y such that Tp{x) — \Tp{r)\. Assuming 
that z < rtv (the opposite case is analogous), we have il){x,r) <OinO<a;<z and 
'4'{x, r) > in in z < a: < r. Bearing in mind that there is a TV of at rtv and 
assuming that 0' passes from postive to negative (the opposite case is analogous), the 
curves of tp^ "0 and 0Vi ©V are qualitatively analogous to those shown in figure p] 
Hence, the proof is the same as that of Proposition 7 with the restriction that r > rtv 
(see Appendix D). 
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Note . While the conditions for C > are similar to those of domains of elliptic 
models with topology M^ examined in Propositions 4-7, these conditions are far more 
restrictive for the topology S'^. This is so because the parameters y and z change from 
one ■^T(i) to the next and are strongly constrained. While we have r/rtv > 1, this 
ratio cannot take arbitrary large values because there is no asymptotic radial range. 
Since regularity conditions provide extra constraints on the profiles of J-' and J-' (see 
Appendix A) , it is practically impossible to make general statements applicable to all 
profiles and all ^T(i). The verification of C > for elliptic models with closed topology 
needs to be done in a case by case basis, though it is likely that regular configurations 
should exist for which Q> holds, at least for some ^T{t). 

13. Elliptic regions in an expanding background. 

The most interesting configurations that can be constructed with LTB models are 
perhaps those describing an inner elliptic region possibly undergoing local collapse 
smoothly immersed in an expanding (parabolic or hyperbolic) background. 



'FJt.x) 



1.0 




jr = rb 



Figure 7. Profiles of J^ and jFp in an elliptic region surrounded by a 
hyperbolic exterior. The figure displays the radial profile of JF and those of 
the time depen dent jFp for various hypersurfaees ^T{t). Since JF = JF(x), the 
profile given by ( |80[ ) is the same for all '^T{t). While the profile of jFp is different 
for differen t ^T (f), a ll sha re common features that emerge from Lemma 3 and 
properties | |58a| and l |58b| l that are valid for all •^T(t). All curves J^p intersect J^ 
when .7-"' = in a value x = y < r^, (white circle to the left), and for all curves 
there is a value x = z > r]^ for which J^p = 1 (white circle to the right). Different 
types of domains r][r] in a elliptic/hyperbolic configuration are displayed in figure 

m 
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Figure 8. Domains in an elliptic region enclosed by a hyperbolic 
exterior. The profile of T{x) is shown (thick curve) next to that of a Tp in 
a fiducial ''^T{t). Assuming absence of TV's of ©, three types of domains ri{r\ 
can be constructed, depending on the location of the domain boundary r (white 
circle). It is evident that domains of the type I comply with \72\ and so are 
identical to elliptic domains without TV's (Proposition 4), those of the type III 
behave as hyperbolic domains without TV's (Proposition 2), while domains of 
type II are identical to elliptic domains with a TV of F (Proposition 6). We can 
then identify three zones: an internal elliptic one (dark gray) and a hyperbolic 
external one (very light grey), with an intermediary transition zone. Notice that 
the boundary of the elliptic region a; = r^ lies in the intermediary zone. If there 
is a TV of 0, we have the same situation, but then, depending on the location 
of the TV, either one of the domains I, 11 and III would be equivalent to one of 
those of Propositions 5, 7 and 3. 
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13.1. Hyperbolic exteriors. 

An elliptic region surrounded by a hyperbolic exterior can be constructed in a single 
coordinate patch (without "matching" as in [2U|53]) by choosing !F as 

{< 1, <x <rb 
= 1, x = rk (80) 

> 1, X > rb 

where rb marks the comoving boundary between the regions. We assume that there 
are no TV's of O and R, though J^(0) = 1, and so JF must necessarily have a TV in 
order to reach T{rh) — 1 (besides the fact that !F and !F' must comply in the elliptic 



region with regularity constraint (1.3 1, see Appendix A). We have a similar behavior 
of J-' and J-'p as in elliptic domains with a TV of J^, except that now the hyperbolic 
background allows for domains in which J- > 1 for a; > rb. The profiles of J^{x) and 
Tp{t,x) for various ^T{t) are shown in figure^ We prove now the following 

Proposition 9: In all regular LTB configurations made by a hyperbolic region 
surrounding an elliptic region containing a center, without a TV of R but 
with possible TV's of G, there are always domains ri[r] with r > a for some 
a > for which C > holds. 

Proof. 



From Lemma 3 and (58a)~(58&l, which are valid for all ^'?"(i), the profiles of Tp in the 
elliptic/hyperbolic configuration are such that in every ^T{t) we have T' = when 
T = Tp, and bo Tp < T and T'p > as a; increases asymptotically (see figure [tI 



We examine the sign of (^ in ( 55 1 for Tp whose profile is as shown in figure It] We 
consider first the case without a TV of Q. Given the properties of the curves of Tp in 
all ^Tif) , there will necessarily exist values y < r^, and z > rb of the radial coordinate 
(see figure It]), such that T{y) = J-p{y) and Tp{z) = 1 hold. It is evident that the 



relation between Tp{x) and Tp{r) in (55 1 for a mixed elliptic/hyperbolic configuration 
strongly depends on the choice of the domain boundary. We use figure [8] to illustrate 
how, given the curves of figure p] the sign of (p depends on the selected ri[r]. 

Following figure [SJ it is evident that the elliptic/hyperbolic configuration allows 



for three types of possible domains 77 [r] and behaviors of </? in (551 along any ^T{t): 

I, < r < y : 

Tp{r)<Tp{x) Va;G?7[r], ^ < 0, C<0 (81a) 

II, y < r < z : 3 yi £ ri[r] such that : 

Tpir)<Tp{x) 0<x<yi, ip<0, (81&) 

Tpir) > Tp{x) yi<x<r, (^ > 0, (81c) 

III, r > z : 

Tp{r)>Tp{x) Va;e77[r] </? > 0, C > 0, (81d) 

Hence, it is evident that C > holds for any domain ri[r] with r > z. Of course, z will 
be different for different ^T{t), but given the fact that for all ^T(i) we have T' > Q 
for X > y, & value z > r\j fulfilling the desired result always exists in every ■^T(t). 

The presence of a TV of Q does not make a significant effect. If this TV occurs in 
domains along the inner or intermediary regions I and II: < x < z, then its effect is 
the same as that discussed in proposition 7 and depicted in the bottom panel of figure 
I5] In this case, C > holds (at least) in the more external part of rj[r] and has no 
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effect on the hyperbolic region III, since would already be monotonous in III, and so 
Proposition 2 would apply. If the TV of Q occurs in the external zone III x > z, then 
we have exactly the same situation as in Proposition 3, and thus, irrespectively of the 
profile of 0, the result of this Proposition would apply: there always exist domains 
ri[r] with r > zi > z for which C > holds. As commented in the proof of Proposition 
3, the presence of several TV's of Q would not alter the result as long as there is a 
clear asymptotic behavior of Q' . 

Note . 

The main implication of Proposition 9 is that the determinant factor for C > is 
the existence of a regular hyperbolic asymptotic exterior, with local features (TV's 
of Q or the enclosed elliptic region) playing a very minor role. In other words, given 
the non-local nature of back-reaction, local features can always be "coarse grained" 
when averaging domains are sufficiently large. As shown in figure [8] all domains of 
the form I are identical to elliptic domains without TV's, while all domains of the 
form III can be treated (as far as conditions for back-reaction are concerned) just like 
hyperbolic domains with or without TV's (of O). Intermediary domains of the form 
II are transitional, being practically identical to elliptic domains with a TV of J-" and 

e. 

Notice that Proposition 9 is valid even if the internal elliptic region undergoes crit- 
ical collapse conditions. If dust layers reach the collapsing singularity, then the '^T(i) 
after some t will no longer be fully regular, becoming singular at the coordinates 
marking the singularity. The range r][r] would be necessarily restricted. However, as 
explained in Appendix A, the function Tp remains regular and this development sim- 
ply requires one to treat the involved integrals as improper integrals, so that 0j,, O^ 
and other scalars are regular at all points save at the singularity locus. 

Note: Spatial curvature in mixed configurations. 

It is interesting to examine the behavior of the spatial curvature ^TZ in "mixed" 
elliptic/hyperbolic regions of this type. Evidently, for positive curvature to become 
negative at a fixed comoving boundary, we must necessarily have a sufficiently large 



^TZ' < 0, which means that the critical value (74 1 is necessarily reached, as in elliptic 
domains with a TV of ^ (Propositions 6 and 7). Since (from (28&I and (61al) a zero 
of 1 — .7-" is a zero of the quasi-local curvature '^TZq and for '^TZ' < we have '^TZ < ^TZ 



then (from (27&)) condition .F(rb) = 1 implies 



9' 



'7^(rb) = -^ / ' ^TZ' V, dx < 0. (82) 

Vq(rb) Jo 

Hence, for values near rb (but < rb) in the elliptic region (for which 1 — JF^ is still 

positive) the local curvature is already negative. This is an example showing that an 

elliptic region (1— JT^ > 0) does not necessarily imply that ^TZ > holds in every point, 

though ^TZ > in every point does imply 1 — J^^ > 0. This situation is illustrated by 

figure [9] 

13.2. Parabolic exteriors. 

If instead of a hyperbolic exterior the elliptic region is surrounded by a parabolic 
exterior, we have 

:Fix) ( <!' 05-<-b (83) 

1 = 1, a; > rb 
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x = 



Figure 9. Spatial curvature in an elliptic region surrounded by a 
hyperbolic exterior. The figure displays t he r adial profile of the spatial 
curvature "^TZ and dual quasi— local ^TZq given by ( 28fe i. The function J^ is diplayed 
for comparison. In order to have a positive ^Ti near the center becoming negative 
at finite comoving coordinate values, there must be a sufficiently l arge negative 
gradient ^TZ' . In these conditions, ^TZ < "^TZq and the critical value | |74[ ) for a TV 
of J^: '^'R, = ^TZq/S can be reached. At this value J^' = and, as x increases, J^ 
increases until it reaches .7^ = 1 at the boundary r = r^, between the elliptic and 
hyperbolic regions. Notice that ^7?. is already negative in areas inside the elliptic 
region where JF < 1 holds. 



It is evident that the replacement of the hyperbolic region by a parabolic one, simply 
removes the external zone (hyperbolic domains of type III) and keeps the intermediary 
zone (elliptic domains with a TV of J^) all the way into the asymptotic range. 

Proposition 10: In all regular LTB configurations made by a parabolic region 
surrounding an elliptic region containing a center, without a TV of R but 
with possible TV's of O, there are always domains rj[r] with r > a for some 
a > for which C > holds. 

Proof . 

Since the behavior of J- and !Fp is qualitatively analogous to that of elliptic domains 
with TV's of T (and possibly O), the proof is afforded by Propositions 6 and 7. 



Note 

It is evident that the elliptic/parabolic configuration is qualitatively analogous to an 
elliptic model with the '^T(i) having topology M.^ and spatial curvature decreasing 
(^TZ' < 0) and ^7?. -^ asymptotically. However, the assumption of the parabolic 
exterior (J-" = 1 for x > r\y) makes it much more artificial and contrived. The most 
interesting result is the fact that back-reaction can be positive in regular regions 
containing parabolic domains only if the latter do not contain a symmetry center. 
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Parabolic domains: J^ = I 


Turning values 


Q and domain restrictions 


Comments 


A>0 


None 
(possibly 6) 


Q = holds for all r][r] 


Proposition 1 


No 


Hyperbolic domains: T > 1 


Turning values 


Q and domain restrictions 


Comments 


A>0 


None 

e 


Q > holds for all r][r] 

Q > holds for r][r] with 
< r < ri or r > r2 with ri < r2 


Proposition 2 
Proposition 3 


Yes 
Yes 


Elliptic domains: < JF < 1 


Turning values 


Q and domain restrictions 


Comments 


A>0 


None 

e 

J^and e 
R 


Q < holds for all r][r] 

Q > holds for 77 [r] 
with < ri < r < r2 

Q > holds for 77 [r] 
with r > a with a > 

Q > holds for 77 [r] with 
< r < ri or r > r2 with ri < r2 

77 [r] with r <rty: as the cases above 
77 [r] with r > rtv: as with TV of J^ 


Proposition 4 
Proposition 5 

Proposition 6 

Proposition 7 

Proposition 8 
Spherical topology 


No 

Extremely 
restricted 

Yes 

Yes 

Extremely 
restricted 


Elliptic regions with a hyperbolic exterior 


Turning values 


Q and domain restrictions 


Comments 


A>0 


J" 
(possibly 6) 


Q > holds for all 77 [r] 
with r > a for a > 


Proposition 9 


Yes 


EUiptic regions with a parabohc exterior 


Turning values 


Q and domain restrictions 


Comments 


A>0 


J" 
(possibly 6) 


as above 


Proposition 10 


Yes 



Table 1. Summary of sufficient conditions for back— reaction in LTB 
models. The domain r}[r] is defined by | |13[ ). Turning values are defined in 
section 8. The parameters ri, r2, rtv and a depend on the profiles of J^ and © 
for each case and are given in the proofs of the Propositions. In elliptic domains 
the regularity conditions strongly constrain the profiles of J^, and so conditions 
stated in the table might fail to hold (see sections 8-12 and Appendices A and 
B). Sufficient conditions for yl > are discussed in the follow up paper (part II). 
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14. Final discussion and conclusion. 

We have provided in this article a fully comprehensive discussion of the back- 
reaction term Q for LTB models in the context of Buchert's averaging formalism. All 
possible regular configurations admitting at least a symmetry center were examined 
in ten propositions that where rigorously proven for parabolic, hyperbolic and elliptic 
domains, regions and models, and for turning values (TV) of 8, J^ and R along 
hypersurfaces ^T{t). A summary of the conditions for Q > is provided in Table 1. 

The form of C for all domains 77 [r] leads to a determined sign of Q for parabolic 
{Q ~ Proposition 1) and for hyperbolic (Q > 0) or elliptic (Q < 0) models without 
TV's (Propositions 2 and 4). In all other configurations, the sufficient conditions for 
the occurrence of Q > depend on the selected domain 77 [r] and on the existence 
of TV's of Q or T. These conditions are not too restrictive for domains and models 
(hyperbolic or elliptic) whose hypersurfaces ^T(i) have topology M'^ and have a full 
radial asymptotic range (Propositions 3, 5-7 and 9-10). For elliptic models whose 
^T(i) have spherical topology (TV of R, Proposition 8), these conditions are much 
more restrictive. 

From the proofs of the Propositions it is evident that the profiles of !F, J-p and 
their radial gradients determine the conditions for C > 0. However, as discussed 



in several parts of the article (for example, see equations (61a|-(61&l in section 8.2, 
as well as sections 10 and 13), these radial profiles are themselves closely related to 
the spatial curvature '^TZ, its average quasi-local function '^TZq and their gradients 
^7?.', ^Ti-q. Hence, statements about configurations favoring Q > can also be given 
as statements on ^TZ in specific domains. Thus, fulfillment of Q > (and possibly 
^cff > 0) is compatible with the following domains 

• ^Tl < for all r ("pure" hyperbolic models of section 10). 

• ^TZ < for large r, even if these domains contain elliptic inner regions with '^TZ > 
(hyperbolic models enclosing an elliptic region, section 13). 

• '^TZ > 0, with sufficiently large gradients ^TZ' < and ^TZ -^ asymptotically 
(elliptic domains with open topology and a TV of T, section 11) 

This connection between the sign of Q and the sign and profile of the spatial curvature 
(specially negative curvature) fits very well with the dynamical importance of this 
curvature due to the presence of large cosmic voids in the pattern of cosmic large scale 
structure [37^. 

However, JF is an invariant scalar in LTB models [5S] and can also be associated 
with a covariant characterization (valid for spherical symmetry) of "binding energy" , 
which has been proposed in [251 EH HH US El] by means of the comparison between 
the quasi-local or "effective" mass-energy function Mq (the function AI in ([3])-(|4])) 
and the proper mass-energy Mp function: 

2Mp{r) = ^ pdVp = K pT-^R^R'dx, (84a) 



(84&) 



IV[r\ 

2Mq{r)^2M{r)^ ^^ I pTdVp^nj pR^R'dx, 

C Jvlr] Jq 

where k — SttG/c^. A binding energy integral functional can be defined as 
B[r] = 2[Mq-Mp] = ^ f p{T-l)dVp = K r p ^^'^l^^^' dx 

C Jv[r] Jo J' 

= ^((^-l)p>V, (85) 
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which for every domain is a conserved quantity along the 4- velocity flow: B[r] — 0. 
Bearing in mind that the radial profiles of J-' and J-'p are qualitatively analogous and 
that B[r] in (85 1 is not a local but an integral quantity associated to a whole domain 



r][r], all statements about a given asymptotic behavior of J-'p can be given as statements 
on the behavior and profiles of this binding energy in different domains. Therefore, 
fulfillment of condition Q > can be associated with domains in which 

• /B[r] > for all r ("pure" hyperbolic models of section 10). 

• B[r] > for large r, even if these domains contain elliptic inner regions with B[r] < 
(hyperbolic models enclosing an elliptic region, section 13). 

• B[r] < 0, where the gradient B' < becomes less negative for large r and B ^ 
asymptotically (elliptic domains with open topology and a TV of JF, section 11) 

This connection between binding energy and back-reaction has been highlighted by 
Wiltshire [2S|- 

It is important to notice that, regardless of the TV's of Q or T, or changes of 
signs of spatial curvature or binding energy and other local complexities of the radial 
profiles, the conditions for positive C (sufficient for Q > 0) are basically determined 
by the asymptotic radial behavior of the incumbent scalars. The dominance of the 
asymptotic behavior is particularly evident in the case of the elliptic region with a 
hyperbolic exterior discussed in section 13 and illustrated by figures [7] and [5] a local 
feature (like the elliptic region or a TV of O) has no effect on the occurrence of C > 0, 
as long as we select domains that lie sufficiently far away from the center. This is 
so, not only in the external hyperbolic region, but even in the transition zone still 
inside the elliptic region (see figure Isl . The fact that the asymptotic behavior is the 
key factor is expected, as back-reaction is a non-local effect, and so local features can 
always be "coarse grained" , even if they exhibit critical conditions like a collapsing 
singularity in the elliptic region. 

The Propositions that we have proved provide information on the existence of a 
positive Q, but not on how large this term can be. It is possible to argue qualitatively 
that, while back-reaction can be positive in elliptic and hyperbolic domains, the latter 
are likely to yield a larger value. This follows from comparing the integrands of <i> and 
^ in (52|-(56l for each type of domain. Assuming the gradient Q' to be of similar 



magnitude in each case, the function ip in (55) is likely to be a small quantity in most 
of the asymptotic radial range of an elliptic domain compatible with C > 0. If we look 
at the profile of JFp displayed by figure 111 as r/y ::^ 1 , we have JF ^ 1 but JF < 1 
for all X. Since J-p < T, then J^p —^ 1 (asymptotically small curvature). But then, 
!Fp{x) w J^pir) w 1, which makes ip in (55 1 close to zero for a long asymptotic radial 
range (notice also that JF and J^' are strongly constrained by ( |1.3| ). In comparison, 
J^p is not restricted to remain below J^p = 1 in the asymptotic range of hyperbolic 
domains or regions (see figures n\ and p| , and so this function can reach much larger 
values than in elliptic regions, which makes the integrand of <i> much less restricted. 

The magnitude of the back-reaction term Q is also related to the magnitude of 
the radial gradients of various scalars. Using (19&I, (21al, (26b), (27a), (49) we can 
write C as 



C = 



R 



el 






-e' 



R 
3R' 



e: 



J^ 



p f 



-el 



e: 



J" 



P -r- 



e: 



(86) 



However, the gradients G' and G' can be related by means of (28c) and (32) to 
the gradients of pp, (p), '^TZq and { TZ), while gradients of ^TZq relate to gradients of 
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J^, then it is evident that C will be larger in radial ranges where these gradients are 
significant {i.e. in regions where inhomogeneity is not negligible). As a consequence of 



(86 1 and the qualitative arguments in the previous paragraph, it seems that hyperbolic 



models that allow for larger gradients of J^ and J-'p (in a wider radial range) are much 
more likely to provide a back-reaction Q that is sufficiently large for yl > in ( 37 1 



to hold for density profiles that are not too restrictive, though even if the resulting 
C > could be small in elliptic models an effective acceleration could still occur if p 
is sufficiently small. 



Since 6p (or (6)) and 6^ in (|28c|) and ([32|) play the role of "kinetic energy" 



6p and 6^ given by (|86|) provides 



terms, the relation between C and the gradients 

theoretical and contextual support to Wiltshire's interpretation of back-reaction [26] . 
Given the LTB configurations which we have proven to be compatible with Q > 0, we 
examine in a follow up paper (part II) the magnitude and radial variation of the back- 
reaction term, as well as boundary conditions that are needed to find the conditions 
for the fulfillment of ( 37 1 and/or ([40|) , which lead to the existence of a positive effective 



acceleration that could mimic the effects of dark energy. Another line for future work 
is to extend the results presented here and in part II to the more general Szekeres 
models [121 ISH]- These tasks are being undertaken in separate continuing articles. 



Appendix A. Regularity of LTB models. 

Standard regularity conditions of LTB models are essential for validity of the results 
presented in this artice. These conditions have been extensively discussed in the 
literature [111 [201 IHJ US! HSl HI] • We provide here a brief summary. 



Symmetry centers 

We have only considered LTB models having (at least) one symmetry center, which 
is a regular timelike worldline corresponding to a fixed point of the rotation group 
SO (3). This is a sufficient condition for integrals in (151 and (23 1 to be finite in a 
domain ri[r] defined by ( 13 1 . The symmetry center can be marked as r = (and r = r^ 
if there is a second one). The following conditions hold: R{t,0) = R{t,0) = 0, while 
.F(O) = 1 and M(0) = M'(0) = T'{0) = 0, but i?' ^ 1 as r -> 0. Also, for any scalar 



A and its dual functions Ap and Aq obtained by means of (15 1 and (23) we have 

Ait,0) = {A)[0] = Ap{t,0) = Agit,0), (1.1a) 

A'{t,0) = {AYiO] = A'p{t,0) = A'q{t,0) = 0, (1.16) 

From ( |4T| , (50 1 and the equations above, it is evident that back-reaction vanishes at 
the center: C(0) = C'(0) = 0. Also, we have from (|42]): A(0,r) = -(k/2)p(0) < 0, 
hence there is a nonzero effective deceleration in all domains sufficiently closed to a 
center in all LTB models for which p{0) > 0. Notice that a central singularity is also 
associated with _R = 0, but its coordinate locus is not a comoving worldline. 



Regularity at a TV of R and shell crossing singularities. 

There is a "shell-crossing" singularity if i?'(i, r) = occurs for coordinates that are not 
of the type r =constant, thus an important regularity condition is then [20l [22ll23ll24] 

i?'(t,r)>0, (1.2) 
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holds with B' = occuring (regularly) only in models whose '^T{t) have spherical 
topology at a TV of R (see Appendix B) . 

Necessary and sufficient conditions to avoid shell crossing singularities can be 
given in full analytic form. This is done in [501 HH US] in terms of the free functions 
associated with ([3| and using the function E related to ^ by JF = [1 + £']^/^, while in 
j24] it is done on initial conditions specified on a fiducial ^T{ti), and using K defined 
by^=[l-X]i/2. 

In hyperbolic domains and models the conditions to avoid shell crossings 
necessarily require M' and JF' to have the same sign as R' , which prevents the existence 
of TV's of JF (see equation (22) of [33]). However, these regularity conditions allow 
for a TV of JF in elliptic models, though the gradients !F' are strongly constrained, as 
shown by equation (23) of [23] and equation (76) of [24|- It is useful to rewrite the 
latter equation in terms of J^ and its gradients, M and its gradients and Ri = R(ti, r) 
where ti marks an arbitrary initial ^T: 



27rR 



jtI™) 



^(k) 



> P^Q^ 



Ava) 



Ak) 



s-(m) c(fc) 



>o, 



where 

.„) _ M'/M 



1, 



TT' 



P, = ± 



SR'jRi 

[l-(l-i;)2]i/2 



U'^^ = 

2 ' {1-T^)R'JR, 



-1, 



Q, = arccos(l - u) - [1 - (1 - vf]^/^, 



(1 - T^) R 

V = 

M 



'Tin 



(1.3) 

(1.4a) 

(1.4&) 
(1.4c) 



Thus, avoidance of shell crossings places much stronger constraints on the radial profile 
of J^ in elliptic models than in hyperbolic or parabolic models. This fact is important 
for the discussion of sections 10 and 11. 



Range restrictions due to a collapsing singularity. 



So far we have assumed that the integration range ry[r] defined by ( 13 1 is fully regular. 
However, it is a well known fact that a collapsing singularity arises in elliptic LTB 
models when < 0, and the coordinate locus of this singularity is (in general) not 
simultaneous (i.e. not marked by a constant t or single ^T{t)). In general, this 
collapsing singularity is marked by a curve [t(rcou), ?'coii] in the (i, x) coordinate plane, 
where R{t{rco\\),rco\\) = and curvature scalars diverge (see [35]). Hence, in any 
collapsing elliptic region the hypersurfaces ^T(t) for t > i(rcou) are only regular for 
the semi open subset 7y[r] = {x | TcoU < a; < r} C ri[r]. However, the existence of 
this singularity has no consequence in the definition of (A) [r] or Ap and Ag functions 
because the involved integrals can be treated simply as standard improper integrals. 
We define at each ^T{t) the incumbent integrals with their lower integration limit as 
y — ''coll + £; for an arbitrarily small e > 0, and then obtain the limit as e — » 0. Off 
course, since & —^ — c» in this limit, {&) or &p or &g might diverge as well, but the 
functions are well defined and behaved in the range fj[r]. Regarding Tp — Vq/Vp, the 
integrals involved it its definition are regular as a; — > rcoii. Therefore, the restrictions 
mentioned above simply prevent J^p from taking values x < rcoii, and so the results 
presented in sections 7-9 for any configuration that involves collapsing layers can be 
trivially extended to include hypersurfaces '^T(t) for t > t(rco\\). 
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Appendix B. Radial dependence and proper length. 

While spherical symmetry effectively reduces the spatial variation of scalars on the 
^T(i) to a one-dimensional radial dependence, there is no inherent covariant meaning 
in the radial coordinate (in fact, the metrics (fTl) and ([9| are invariant under arbitrary 
re-scalings r = r(f)). However, radial rays are totally geodesic (they are spacelike 
geodesies of Q), so it is possible and desirable to relate radial dependence at each 
^T(i) to dependence on the natural (and covariant) parametrization of the rays in 
terms of their affine parameter or proper length: 

e(r) = / y^dx = / ^dx, (2.1) 

Jo Jo -^ 

so that ^(0) = 0, and we have adopted the notation L = J^~q , so that (unless 
specified otherwise) all functions inside an integral depend on the dummy variable x 
(and as stated before, i is a fixed parameter). 

Since f is a proper length, we must necessarily have ^(r) > and f'(r) > for 
r > 0, so that r2 > ri implies ^(r2) > S.{ri) and the converse is also true. These 
properties define through pTTj) a well behaved radial coordinate as that complying at 
every ^T(t) with the following regularity condition 

sign (i?') = sign J^, 3 TV of i? ^ 3 zero of J^, (2.2) 

where "TV of i?" stands for a regular "turning value" of R, or a zero of the gradient 
R' under regular conditions (see section 8) . [^ 

The relation between profiles of scalars with respect to r and of ^ at each ^T is 
given by 

,, dA ^, dA R' 

A' -- f -- 

dR dRR^ 



which (as long as (2.2 1 holds) implies that 

sign {A') = sign {dA/dC) , 3 TV of ^(r) 

sign (i?') = sign {dR/dO, 3 TV of i?(r) 

so that radial dependence is qualitatively analogous to proper length dependence at 
each ^T(i): monotonous radial profiles A{r) correspond to a monotonous profiles A{£^) 
and a TV of A corresponds to a zero of dA/dS,- Notice that (|2.36l allows us to relate 





(2.3a) 


OR 


(2.36) 


^ 3 TVof A(0, 


(2.4«) 


^ 3 TVofi?(0, 


(2.4&) 



J^ to initial value Ri — R{ti,r) at at a fiducial hypersurface^^. Since this relation 
between T and dR/d^ is valid in all ^T(t), a zero of J^ (characteristic of S^ topology) 
or a monotonous profile of J- (characteristic of M'^ topology) will be common to all 
^T{t). 

Appendix C. The fluid— flow evolution equations. 

LTB models are usually examined by means of the solutions of the Friedman-like 
equation (pi). Since these models (as all spherically symmetric spacetimes) can be 

^ This definition obviously excludes shell crossing singularities discussed in the previous Appendix 



for which R' = occurs in violation of ( 2.2 I 
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completely described by covariant scalars |30j . an alternative approach to study their 
dynamics is through the evolution equations for the scalars 

^ = {p, e, s, £, 37^}, (3.1) 

which follow from equations Q and (|5|-(I7| and provide a complete set to characterize 
the LTB models. This is the "fluid flow" or covariant "1+3" framework of Ehlers, Ellis, 
Bruni, Dunsbury and van Ellst [371 UHl IMl 130] and leads to the evolution equations 

92 



6 



- p - 611^ 



p = 


-pQ, 




t = 


--...- 


-£, 


£ = 


-fpS-35 


fe 



together with the spacelike constraints 



3i; 



R' 



0: 



6 



p 



S' + 3£ 



R' 



R 



0, 



6^ 

yy ' "~ R 

and the "Hamiltonian" constraint (analogous to the Friedman equation) 

2 ,^ 3t 



= 3^ 



6 



(3.2a) 
(3.26) 
(3.2c) 

(3.2rf) 



(3.3) 



(3.4) 



which is an integral of the Raychaudhuri equation (3.2a I 



Buchert's equations (30)"(35l follow directly by applying the proper average 



functional (15) to both sides of equations (3.2a), (3.2&I and (3.4) (the Raychaudhuri 



and energy balance equations and Hamiltonian constraint), and using the time 



derivative rule (20), as well as the formulae for the covariance and variance momenta 



(18&I and (18c I. The integrability condition 

Q + 2(e) Q + lie) Cn) + (^n) • ^ o, 



(3.5) 



can be proven to be compatible with the fluid flow evolution equations (3.2a|-(3.2dl 



and the Hamiltonian constraint (3.4). The proof follows after a very long algebraic 



and (18c) 



manipulation expressing E and £ as in (29a|-(296), together with applying (20 1, (18&I 



Appendix D. Proofs of Propositions 7 and 8 

Appendix D.l. Proposition 7. 

Consider an elliptic domain with a TV of J^ and a TV of e (but not of R). The 
profiles of e'(f and O'-0 for this configuration are those displayed in the bottom panel 
of figure pi) , with x — y marking the location of the sign change of f and x = rtv 
marking the TV of 0. This leads to 

$ = - / e'\Lp\Ax+ / e'ipdx- / \e'\Lpdx, (4.1a) 



* = 



e'V'dx- / le'lV'da;., 

Jrtv 



(4.16) 
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We introduce now the real positive numbers {a, (3, 7, J, e} by 
rv /-'tv 

0</ &\ip\dx<ay, 0< Q'<fdx<f3{rty-y), (4.2a) 

Jo Jy 

0</ |e'|(/3da;<7(r-rtv), 0</ O' ^/j dx < 5rt^ , (4.26) 

Jrtv "'0 

0</ |e'|?/'dx<e(r-rtv). (4.2c) 

Condition C = $$ > implies 

< mm < 1 1 , 1 H — > (4.3) 



rtv L 7 ''tv 7 e 

or 

— >maxn -^^ + -,1+-L (4.4 

Ttv L 7 ^tv 7 ej 

which provides the values of the real numbers r\ and r2. 

Appendix D.2. Proposition 8. 

The profiles of 0'(p and Q''tp associated with Proposition 8 are qualitatively analogous 
with those of figure |5] and so the proof above applies also to that case. The only differ- 
ence is that the role of y above corresponds to the value z such that J-p{z) = |J^p(r)|, 
which marks a zero of tp where this function passes from negative to positive (see 



(79a|-(796) in section 12). However, since r/rtv cannot take arbitrary large values, 



the fulfillment of condition (4.4) is much more restrictive 



Appendix E. Proof that a TV of R must be a common TV of Q. 

A TV of R {i.e. R'{rtv) = for rtv G vi^]) occurs in regular conditions only in elliptic 
models whose ■^T(i) are homeomorphic to S'^. We prove in this section that 

i?'(r-tv) = => e'(rtv)-0, (5.1) 

and show (as part of the proof) that this implication also holds for all quasi-local 
scalars Aq obtained from (|23|, as well as for the scalars A = p, ^TZ, M, but not for 



^, J-p and scalars Ap that follow from (15). Notice that the converse statement of 



(|5.1| is false: there can be a TV of 6 with R' > 0. 

From the regularity condition (2.2), a well defined proper radial length (|2.1| 



requires the ratio R' /T to be continuous (at least C") at r = rtv Since i?'(rtv) — 
^{rtv) = 0, with R"{rtv) < 0, J-"'(rtv) < 0, then from I'Hopital rule we have at the 
limit r — > rtv 

which implies that rtv is also a TV of the quasi-local volume Vq (because V'q — 
A-kR^R'), but not of the proper volume Vp (because Vp = AttR^R' /T). By considering 



(5.2), together with (196), (21al and (|27a| 



A'q^^-^[A~Aq]^ A;=^[A-Ap], (5.3) 
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i?'(r-tv) = => A;(rtv) - 0, 
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(5.4) 



but, because of (5.2 1, rtv is (in general) not a TV of the p-scalars Ap. 

To prove (5.1 1 we assume that rtv is a TV of R and examine the behavior of Q 



and 0g around r = rtv + efor |e| ^ 1. Considering (from (5.4|) that a TV of i? implies 
8' (rtv) = 0, we have at leading orders: 

e,(rtv + e) « e,(rtv) + ^e;'(rtv)e', i?(rtv + e) « i?(rtv) + ^i?"(rtv)e', (5.5) 



Since, in general, O(rtv) ^ 0, by applying (5.3 1 to Q and O, and at leading orders we 
get 



e(rtv + e) - e(rtv) = e,(rtv + e) - e,(rtv) + 



3R' 



nv+c 



' %R 

'SR' 



(5.6) 



Using (5.5 1 to evaluate this expression at leading orders and dividing both sides by e 
we get 

e(rtv + e) - e(rtv) _ eg(rtv + e) - e,(rtv) _^ 1 

6 



+ -e;'(rtv)e. (5.7) 



e e 

Taking in both sides the limit as e ^ we get G'(rtv) — ©^(''tv) = 0, which is the 
desired result. This proof can be extended to other scalars (except JF which necessarily 
vanishes at rtv). 
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